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Collective molecule formation in a degenerate Fermi gas via a Feshbach Resonance 
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We model collisionless collective conversion of a degenerate Fermi gas into bosonic molecules via 
a Feshbach resonance, treating the bosonic molecules as a classical field and seeding the pairing 
amplitudes with random phases. A dynamical instability of the Fermi sea against association into 
molecules initiates the conversion. The model qualitatively reproduces several experimental obser- 
vations [Regal et al., Nature 424, 47 (2003)]. We predict that the initial temperature of the Fermi 
gas sets the limit for the efficiency of atom-molecule conversion. 
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PACS numbers: 03.75.Ss, 03.75.Kk, 03.65.Sq, 05.30.Fk 

The idea [l| that an adiabatic sweep across an atom- 
molecule resonance can transform an atomic condensate 
into a molecular condensate has recently been ported to 
experiments on degenerate Fermi gases. By sweeping 
a magnetic field across a Feshbach resonance, at least 
part of the atoms have been demonstrably converted 
into molecules J2J, y] ■ Magnetoassociation of atoms into 
molecules via a Feshbach resonance is also the key to ex- 
periments in which formation of a molecular condensate 
out of a degenerate Fermi gas has been observed J4J . 

To date, most experiments on magnetoassociation of 
fermionic atoms into molecules have been done in the 
collision-dominated regime. Collisions induce thermal 
equilibrium, and statistical mechanics, or indeed thermo- 
dynamics, seems to be the appropriate theoretical frame- 
work p (see also Cubizolles et al. y|)- An obvious excep- 
tion is the adiabatic-sweep experiments of Ref. [2| on fast 
enough time scales that particle collisions are not a major 
factor [2| • These experiments are the domain of the time- 
dependent Schrodinger equation. Time-dependent asso- 
ciation of a Bose-Einstein condensate of atoms into a con- 
densate of molecules has attracted much interest |l|, |7| . 
Treating the condensates as classical fields, as opposed 
to quantum fields, guides and simplifies the analysis of 
Bose systems. In contrast, even the zero-temperature 
Fermi sea of atoms presumably cannot be represented as 
a classical field, a "macroscopic wave function" . This is 
the technical dilemma we set out to tackle. 

Here we develop a collisionless model for magnetoas- 
sociation of a two-component Fermi gas into bosonic 
molecules, treating the boson field classically. In this set- 
ting atom-molecule conversion builds up from a dynam- 
ical instability. We report on comparisons with experi- 
ments [2j , and make the testable prediction that tempera- 
ture limits the conversion efficiency in an adiabatic sweep 
of the magnetic field across the Feshbach resonance. 

Consider a free two-component Fermi gas (spin-up and 
spin-down) with the annihilation operators c k -[ and c k j 



for states with momentum ftk, and the corresponding 
Bose gas of diatomic molecules with annihilation opera- 
tors 6k- Absent collisions, the Hamiltonian reads 
tj ^ r 

^ ^ [e k (4 T c kT + cl iCki ) + (6+ ie k )6 k 6 k 
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Here Se k = Ue^ = h 2 k 2 /2m is the kinetic energy for an 
atom with mass m, hS is the atom-molecule energy dif- 
ference that is adjusted by varying the magnetic field, 
and K kk ' are matrix elements for combining two atoms 
into a molecule. For s-wave processes, K kk ' are func- 
tions of the relative kinetic energy of an atom pair. The 
Hamiltonian Q conserves the invariant particle number 
N = X)k( c kT Ck T + c ki Ck l + 26J c 6 k ). Given the quantiza- 
tion volume V, the invariant density is p = N/V. 

In the spirit of classical field theory, the main assump- 
tion of our model is that the boson operators in the 
Heisenberg equations of motion are declared to be com- 
plex numbers. To facilitate the numerics, we furthermore 
only keep the molecular mode with k = 0. We also as- 
sume that initially the occupation numbers of the spin-up 
and spin-down fermions are the same, and that the sam- 
ple is rotationally invariant. The expectation values of 
the relevant operators then depend only on the energy of 
the state k, he = h 2 k 2 /2m. We scale the c-number molec- 
ular amplitude as j3 = y/2/N bo, define the fermion occu- 
pation numbers P(e) = (c^c^) — (c k ,c k |), the pairing 
or "anomalous" amplitudes C(e) = (c_ k |C k |), and find 
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2eC(e) + ^nf(e)P[l-2P(e)], (2) 

-^Q/(e)[/3C*(e)-/3*C(e)], (3) 

W + 7r4^ I deV~ef(e)C(e) . (4) 
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Here hep = h 2 (3n 2 p) 2 ' 3 /2m is the usual Fermi energy. 



The energy-dependent atom-molecule coupling is «(e) = 
k(0)/(c) with /(0) = 1, and the Rabi-like frequency is 
fl = >/Nk(0). As per Javanainen and Mackie [l|, k(0) oc 
1/y/V, so that fi ex y/p. The integral arises from the 
continuum limit of the sum over k, and y/e is a phase 
space factor responsible for the Wigner threshold law for 
the dissociation rate of molecules into atoms. 

The problem with Eqs. (J2EJ) is that (3(t) = C(e, t) = 
is always a solution. However, this solution may be un- 
stable. To illustrate, we carry out the linear stability 
analysis of Eqs. 121 II) around the trivial solution for given 
occupation numbers P(e). In the usual single-pole ap- 
proximation that becomes increasingly accurate in the 
limit O — > 0, and for 6 > 0, the Fourier transform of the 
small deviation from j3 — C(e) = has the pole 
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If the fermion occupation number satisfies P(e) 
some energy he, for a suitable detuning S the evolution 
frequency has a positive imaginary part. The implication 
is that a small perturbation from the stationary state 
(3 = C(e) =0 grows exponentially. 

If dissociation of an isolated molecule into two atoms is 
energetically possible, it will invariably happen because 
the state space for two atoms is much larger than for 
a molecule. On the other hand, a filled Fermi sea of 
atoms may block dissociation. The state space of allowed 
molecular states is then the one that is larger, and the 
atoms are prone to spontaneous magnetoassociation into 
molecules. This is the nature of the instability. 

The Fermi sea is thermodynamically unstable against 
formation of Cooper pairs [8(, and resonance superflu- 
ids |9j inherit an analog of this trait of BCS supercon- 
ductors. Nonetheless, suggestive as the similarity may 
be, the BCS instability is different from the present one. 
Thermodynamical instability occurs because pairing low- 
ers the energy, and so coupling to a reservoir with a low 
enough temperature leads to pairing. The hallmark of 
dynamical instability is that a small perturbation grows 
exponentially in time, environment or no environment. 

Quantum fluctuations could trigger spontaneous mag- 
netoassociation, but they are absent in our model. While 
much is known about modeling of quantum fluctua- 
tions classically in boson system |l0j . no correspond- 
ing general methods exist for fermions. We resort to 
the following heuristic device. Instead of starting out 
a calculation with zero anomalous amplitudes, we ini- 
tially seed them with random numbers having zero aver- 
age; specifically, nonzero numbers with a random phase, 
(c- k |C kT ) = e^(4 T c kT ) 1 /2( c j iiCki )i/2. We then intc _ 

grate Eqs. J5J- J3J, whereupon the initial instability and 
the subsequent dynamics run their courses. We do the 
calculations for many choices of the random phases </>k, 
and average the results. This procedure correctly repro- 
duces the initial evolution ex t 2 of the expectation value of 
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FIG. 1: Energy of the bound state of the dressed molecule huj 
as a function of the magnetic field B from the experiments 
of Ref. J2§ (filled circles), and from our calculations using the 
best-fit parameters M, S, and A/x (solid line). 



the number of molecules in the full quantum model. The 
approximation lies in using it at all times. As a technical 
detail, in our numerical calculations we discretize P(e) 
and C(e) at equidistant points of e separated by Ae, and 
resort to the analogous process. 

Next we estimate the parameters of the model. First, 
atom-molecule conversion is unlikely to be the result of 
a contact interaction. Given a nonzero range, there is a 
cutoff in momentum/energy for the matrix element Kkk' ■ 
We crudely assume that the coupling between atoms and 
molecules follows the Wigner threshold law up to the 
point when it abruptly cuts off at some energy KM , and 
correspondingly set /(e) = 9(M — e). Second, we write 
the dependence of the detuning on the magnetic field 
as 6 — A/i(P — B' )/h, where B' is a tentative posi- 
tion of the Feshbach resonance and A/x stands for the 
thus far unknown difference of the magnetic moment be- 
tween a molecule and two atoms. Third, we have an 
atom-molecule coupling strength with the dimension of 
frequency, S, defined by Q = S 1 / 4 S 3 / 4 y / p/m 3 / 4 . 

We ignore the Fermi statistics by setting P(e) = in 
Eqs. (IL'I II) . This yields a linear description of the cou- 
pling between molecules and atom pairs. With A/x > 
and for detunings less than a threshold value So, the re- 
maining set of equations has a stationary solution; the 
Fourier transforms f3(w) and C(e, ui) have a real pole at 
a frequency v(5) such that ui(So) = and u(5) < for 
5 < Sq. What in the absence of the coupling ex $7 were a 
"bare" molecule and a pair of atoms become a "dressed" 
molecule. We interpret faj[8(B)\ as the energy of the 
bound state of the dressed molecule for the given mag- 
netic field, and shift the value of the resonance field from 
B' Q to Bq so that (So = 0. The "renormalized" Bq should 
equal the position of the Feshbach resonance in the limit 
of a dilute gas. 
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FIG. 2: Experimental Uj (filled circles) and simulated (solid 
line) numbers of atoms remaining when the magnetic field 
is swept across the Feshbach resonance toward lower values. 
The simulations use the known experimental parameters and 
parameter values fitted as in Fig. Q 



Finally, we fit the unknown parameters M, A^t, and 
S to best reproduce the experimental binding energy of 
the molecule reported in Fig. 5 of Ref. |2|. The fit min- 
imizes the relative error between the calculated and the 
measured values. The parameters are M = 2ir x 100 kHz, 
Afi = 0.19 fi B , and 5 = 2tt x 580 MHz. 

By setting P(e) = we have ignored the many-body 
shift of the Feshbach resonance. We could include the 
shift by allowing P(e) ^ 0. For our fitted parameters, 
experimental densities, and untouched zero-temperature 
Fermi sea, the shift could exceed one Gauss toward the 
direction of high B. The preparation of the sample in 
the experiments |2| quoted in our Fig. Q] alters P(e) and 
modifies the shift in a manner that is difficult to account 
for self-consistently. Nonetheless, if we insert a full one- 
Gauss shift by hand and fit again, the most relevant pa- 
rameter fl increases by just 60 %. Including or excluding 
the many-body shift of the Feshbach resonance should 
not make a qualitative difference in the values of the fit- 
ted parameters. 

Armed with the parameter values, we next simulate a 
sweep of the magnetic field across the Feshbach resonance 
as in Fig. 1(a) of Ref. |2| by integrating Eqs. (121 H . From 
the experimental maximum density p = 2.1 x 10 13 cm -3 
we have tp — 9.2 x 27rkHz. The initial atomic occupation 
numbers P(e) are set according to the experimental tem- 
perature, ksT/heF = 0.21. The discretization step for 
the atomic occupation numbers and anomalous ampli- 
tudes is Ae = ep/100. We run the magnetic field sweep 
64 times for different random phases of the anomalous 
amplitudes, and average the results. FigureElshows both 
the measured atom numbers (filled circles) and our cal- 
culations (solid line) as a function of the magnetic field 
when it is swept at the rate (40/xs/G) _1 . 
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FIG. 3: Fraction of atoms converted into molecules in a 
sweeps such as in Fig. |21 except that here the sweep rate is 
slower, (400 fis/G) -1 , and the results are plotted as a function 
of the varying initial temperature of the atomic gas. 



Experimentally, half of the atoms are converted into 
molecules at (40/xs/G)" 1 , while our calculations give a 
30 % conversion. There are a number of reasons why 
a full quantitative agreement cannot be expected, e.g., 
keeping only one molecular mode, using the maximum 
density of the atoms instead of making an allowance for 
the density distribution in the trap, simplistic modeling 
of the energy dependence of the atom-molecule coupling 
/(e), and ignoring many-body shifts while fitting the pa- 
rameters. Nonetheless, our model seems to identify the 
correct physics parameters, first and foremost the Rabi- 
like frequency fi, and gives a reasonable estimate for their 
values. 

In the experiments |2j the magnetic field was also swept 
back and forth, whereupon the molecules all dissociate 
back into atoms. With the random initial phases of the 
anomalous amplitudes it is not obvious that our simula- 
tions should reproduce this feat, but they do. 

The puzzling feature of the experiments [2j is that, 
no matter how slow the sweep rate, the conversion effi- 
ciency is limited to about 50%. To investigate, we carry 
out magnetic field sweeps as in Fig. |2 except that the 
sweep rate is such that the detuning as a function of 
time behaves as 5 = — ££1 2 £ with £ = 0.05, corresponding 
to B = (400/is/G) _1 . If fi _1 determines the time scale, 
for £ = 0.05 one would expect adiabaticity. We vary the 
temperature, and plot the final conversion efficiency as 
a function of the temperature. The results are shown in 
Fig. 13 At T = 0, 98% of the atoms are converted into 
molecules. However, by the time the temperature has 
reached the Fermi temperature, kgT — hep, the conver- 
sion efficiency has dropped to 2%. At the typical ex- 
perimental temperatures with ksT ' /hep ~ 0.2 — 0.3, the 
efficiency in fact hovers in the neighborhood of 50%. 

For an increasing temperature, a decreasing number of 
the initial atomic states have a thermal occupation num- 



ber of at least |. The instability responsible for initiating 
atom-molecule conversion then acts on fewer and fewer 
atomic modes. We presume this is the reason for the 
temperature dependence. 

The weakest point in our argument may be the as- 
sumption of a single molecular mode, according to which 
the molecules emerge as a condensate. It would appear 
that, without any seed for the molecular condensate or 
initial phase from the atoms, the atoms are converted 
into a coherent molecular condensate in a time that is 
independent of the size of the sample. This can hardly 
be correct. Inclusion of a multitude of molecular modes 
would probably cure this shortcoming. By construc- 
tion we would still have a classical field representing the 
molecules, but it could consist of patches with uncorre- 
cted phases. A numerical modeling of this situation is in 
principle possible, but much more demanding than our 
present calculations. So far we have made little progress 
in this direction. 

Instead, we present an estimate for the size of the 
patches. The natural velocity parameter in this problem 
is the Fermi velocity, Vf = y/^hep/ra. Suppose the con- 
version takes place in a time At, then atoms in a region 
of size £ ~ vf At would plausibly be able to form a patch 
of molecular condensate. Let us estimate the conversion 
time from Fig.|3 say, as the difference between the times 
when the number of molecules rises from j to | of its 
final value in the dashed-line data, then the size of the 
patch and the characteristic distance between the atoms 
are related by £ ~ 2 / o~ 1 / 3 . Our modeling should be ade- 
quate if it is taken to represent a patch of about 8 3 ~ 10 
atoms. Within each patch we could still resort to the sin- 
gle mode-approximation and use the Bose-enhanced cou- 
pling strength Q that depends on density but not directly 
on the number of atoms. However, the patchy molecu- 
lar condensate would not appear as a condensate in an 
experiment, but would display a momentum distribution 
not narrower than ~ H/£. 

The number of atoms in a zero-temperature Fermi sea 
is the same as the number of occupied states, so for con- 
sistency we should use a step of the order of Ae ~ ejr/10 
in our modeling. We have used Ae = cf/100 in Figs. 
and[3J Now, constants of the order of unity are beyond 
our dimensional- analysis argument, and the step size es- 
timate should be taken with a grain of salt. On the other 
hand, numerically, the step size dependence is logarith- 
mic and weak. Simulations with just one molecular mode 
should therefore be valid semiquantitatively. 

Starting from a microscopic many-particle theory, we 
have modeled collisionless collective conversion of a de- 



generate Fermi gas into bosonic molecules via a Feshbach 
resonance. The key techniques are to treat the molecules 
as a classical field, and to seed pairing amplitudes with 
random phases. The main concept we have unearthed is 
dynamical instability of a Fermi sea against magnetoasso- 
ciation into molecules. The same instability should occur 
in photoassociation. Our model reproduces qualitatively 
all experimental observations |2| where we have tried a 
comparison. Moreover, we have the testable prediction 
that temperature sets the limit for the efficiency of atom- 
molecule conversion. 
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